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We develop a general theory for percolation in directed random networks with arbitrary two 
point correlations and bidirectional edges, that is, edges pointing in both directions simultaneously. 
These two ingredients alter the previously known scenario and open new views and perspectives 
on percolation phenomena. Equations for the percolation threshold and the sizes of the giant 
components are derived in the most general case. We also present simulation results for a particular 
example of uncorrelated network with bidirectional edges confirming the theoretical predictions. 
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I. INTRODUCTION 

A wide class of real systems of many interacting ele- 
ments can be mapped into graphs or networks. Under 
this approach, vertices or nodes of the network represent 
the elements of the system whereas edges or links among 
them stand for interactions between different elements. 
This mapping has triggered a huge number of works and 
a surge of interest in the field of complex networks that 
has lead to a general framework within which to analyze 
their topology as well as the dynamical processes running 
on top of them 0,1113. In many cases, these dynamical 
processes are directly related to functionality and involve 
some kind of transport or traffic flow. Furthermore, the 
very existence of those networks could be naturally ex- 
plained as a direct consequence of the communication 
need among its constituents. The Internet or the World 
Wide Web are clear examples 0. In order to preserve 
functionality, networks characterized by transport pro- 
cesses must be connected, that is, a path must exist be- 
tween any pair of nodes, or, at least, there must exist 
a macroscopic portion of vertices -or giant component- 
able to communicate. In this context, percolation theory 
appears as an indispensable tool to analyze the condi- 
tions under which such connected structures emerge in 
large networks. 

The general theory of percolation phenomena for un- 
correlated undirected random networks was first devel- 
oped by Newman et al. Q after previous results in 
0, . The phase transition at which the giant compo- 
nent forms was well characterized, and the size distri- 
bution of connected finite components below and above 
the critical point were calculated as well. Some further 
refinements were needed in order to approach real nets. 
Hence, correlations between degrees of neig hboring ver- 
tices were taken into account in [l^, UM> ^^'^ gro wing 
networks were treated by Dorogovtsev et al. [13 and 
Krapivsky et al. 

To go further, directness must be taken into consider- 
ation since some of the most interesting real networks 
present asymmetric interactions. Noticeable examples 
are the World Wide Web ' citation networks Ilfil[l7l| , 
email networks gene regulatory networks 'l!^, or 



metabolic networks pOj. Percolation theory for purely 
directed networks was first developed by Newman et al. 
and later by Dorogovtsev et al. jU . In the par- 
ticular case of scale-free degree distributions a number 
of interesting specific results were obtained ■ While 
allowing for general correlations between the incoming 
and the outgoing number of edges of a given vertex, all 
these studies refer to networks with no degree correla- 
tions. Furthermore, the theory is restricted to the class of 
directed networks with no bidirectional edges, although 
this class of edges are ubiquitous to all real directed net- 
works (see 23| and references therein). All these limit 
the applicability of the theory to real networks since bidi- 
rectional edges and degree correlations are common to all 
real directed networks. In this paper we present a general 
theory for percolation in directed random networks with 
general two point correlations and bidirectional links. We 
will show that the presence of bidirectional edges and 
degree correlations modify the picture previously drawn 
6, 7, 21, 22] in a non trivial way, opening new scenarios 
for percolation phenomena. 

The paper is organized as follows. In section II, we 
review concepts, definitions, and the main results previ- 
ously obtained in the analysis of percolation in random 
networks. In section III, we develop the general theory 
for directed networks with bidirectional edges and ar- 
bitrary degree correlations and we show how the theo- 
ries for undirected and purely directed random networks 
stand as particular cases. Section IV is devoted to the un- 
correlated case, which deserves special attention as a null 
or benchmark model. The relative sizes of the giant com- 
ponents are computed and the explicit expression for the 
percolation condition is provided. The well-known crit- 
ical points signalling the phase transition in undirected 
and purely directed uncorrelated networks are recovered 
as limiting cases. A practical application of the formal- 
ism for uncorrelated networks is presented in section V, 
where the transformation from a purely directed network 
to a purely bidirectional one is studied. For power-law 
degree distributions, the transformation is shown to un- 
dergo a nontrivial phase transition. Simulation results 
support this prediction, finding an excellent agreement 
with numerical solutions of the theoretical equations. Fi- 
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nally, we conclude with a brief report of results in section 
VI. 



II. PERCOLATION IN UNCORRELATED 
PURELY DIRECTED NETWORKS 

The topological structure of directed networks is more 
complex than that associated to undirected graphs. The 
edges associated to each node in a directed net are usu- 
ally differentiated into incoming and outgoing. Usually, 
no bidirectional links are considered so that each vertex 
has two coexisting degrees, ki and ko, which sum up to 
the total degree k — km + fco«t- Hence, the degree dis- 
tribution for a directed network is a joint degree distri- 
bution P{ki, ko) of in- and out-degrees, which in general 
may be correlated. The bidirectional edge symmetry of 
undirected networks is thus completely broken in purely 
directed ones, with implications down to the level of per- 
colation properties. The giant connected component in 
undirected graphs becomes internally structured in the 
case of directed networks so that four different types 
of giant components may arise. Whether giant or not, 
these components are characterized as follows (according 
to definitions in [IJ): 

• The weakly connected component, WCC, the per- 
colative cluster in undirected graphs. In the WCC, 
every vertex is reachable from every other, provided 
that the directed nature of the edges is ignored. 

• The strongly connected component, SCC, the set 
of vertices reachable from its every vertex by a di- 
rected path. 

• The in-component, IN, all vertices from which the 
SCC is reachable by a directed path. 

• The out-component, OUT, all vertices which are 
reachable from the SCC by a directed path. 

Notice that, with these definitions, the SCC in included 
in both the IN and OUT components. The percolation 
theory developed by Newman et o/.0,|3| and Dorogovtsev 
et al. for directed graphs with arbitrary degree distri- 
bution and statistically uncorrelated vertices has shown 
that there are two phase transitions: the one at which 
the giant weakly connected component(GWCC) appears, 
and the one at which the other three giant components 
appear simultaneously: the giant in-component(GIN), 
the giant out-component(GOUT) and the giant strongly 
connected component(GSCC) as the intersection of the 
other two. The first phase transition corresponds in fact 
to the standard phase transition in an undirected ran- 
dom graph with arbitrary degree sequence and statisti- 
cally uncorrelated vertices. The condition for thisphase 
transition was first given by MoUoy and Reed @| and 
reads 

^fc(fc-2)P(/c) > 0. (1) 

k 




FIG. 1: Schematic representation of the component structure 
of a directed network. 

When this condition is fulfilled, and although some dis- 
connected finite components may remain, the GWCC 
emerges. It contains a macroscopic portion of the vertices 
in the network capable to communicate to each other re- 
gardless of the orientation of their links. The second 
phase transition is characteristic of directed networks. 
The critical point 

^ fco(fc, -l)P(fc„fco) = 0, (2) 

marks the first simultaneous appearance of the other 
three giant components: the GSCC, the GIN and the 
GOUT, as well as other secondary structures such as 
tubes or tendrils 

The efforts to understand how this landscape is mod- 
ified by the consideration of correlations between de- 
grees of neighboring vertices has been exclusively fo- 
cused on the percolation analysis of undirected networks 
[13 mill- Assortative mixing by degree, observed in 
the vast majority of social real networks, has been found 
to favor percolation in the sense that the giant compo- 
nent appears at lower edge density. On the contrary, 
disassortative correlations, characteristic of technologi- 
cal and biological networks, difficult the formation of the 
giant component even if the second moment of the degree 
distribution diverges. 

III. GENERALIZED PERCOLATION 

In this section we develop a general theory for perco- 
lation in directed random networks with arbitrary two 
point degree correlations and bidirectional edges. We 
term our networks "random" -or Markovian- because, 
apart from purely local properties and two point correla- 
tions being fixed, networks are maximally random. This 
implies that the whole topology is encoded into the de- 
gree distribution P(k) = P(fci, kg, kb) and the transition 
probabilities Pi(k'|k), Po(k'|k), and Pfc(k'|k) measuring 
the probability to reach a vertex of degree k' leaving from 
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a vertex of degree k using an incoming, outgoing, or bidi- 
rectional edge respectively. Notice that we have used the 
notation k = (ki,ko,kb) and that we consider all three 
kind of edges as independent entities. These transition 
probabilities are related through the following degree de- 
tailed balance conditions 

fcbP(k)Pb(k'|k) = fc;P(k')P6(k|k') (3) 

and 

fc<,P(k)F„(k'|k) = fc,:F(k')P»(k|k'). (4) 

These conditions assure that any edge leaving a vertex 
points to another vertex or, in other words, that the net- 
work is closed. Notice that this may not be the case in 
situations where information is incomplete and the out- 
degree of a vertex is known but not the neighbors at the 
end of these edges. 

To analyze the percolation properties for this class of 
networks, it is necessary to calculate the joint distribu- 
tion G{s,s') of the number of vertices (plus itself), s, 
that are reachable from a given vertex and the number 
of vertices (plus itself), s', that can reach that vertex. 
Notice that we can leave a vertex using an outgoing or 
bidirectional edge and we can arrive to a vertex through 
an incoming or bidirectional edge. These sets of vertices 
are called the out- and in-components of a given vertex, 
respectively. Analogously, we can define the marginal 
probabilities Go{s) and Gi{s) for the number of reach- 
able vertices from a given one and the number of vertices 
that can reach it, respectively. These three functions 
contain the information of the sizes of the different giant 
components of the network. Notice that, if the network 
is above the percolation threshold, Go(s) ^ J2s' 
and Gi{s) ^ J2s'^(^'^^)- '^^^ probability to belong 
to a finite component is smaller than one in this situ- 
ation, which implies that the remaining corresponds to 
the probability to belong to an infinite component, that 
is, the giant component. Function Go(s), for instance, 
measures the probability that a given vertex has a finite 
out-component of size s regardless of the size of its in- 
component, which can be finite or not. On the other 
hand, function J2s' ^(■^i only accounts for finite com- 
ponents. Therefore, the relative sizes of the different gi- 
ant components of the network can be written as 

GOC/r= l-^G,(s) , G/iV = l-^Go(s), (5) 



where go{s\]i.) (^^(slk)) is the distribution of the number 
of reachable vertices from a vertex given that we have 



and 

GSGC ^l-Y, Go{s) - G^{s) + J2 (6) 

s s s,s' 

where the GOUT is thought of as the set of vertices with 
an infinite in-component, the GIN as the set of vertices 
with an infinite out-component, and the GSCC as the set 
of vertices with infinite in- and out-components simulta- 
neously. 

In heterogeneous networks, this set of probabilities de- 
pend on the degree of the vertex from where we start the 
count. Therefore, functions G(s,s'), Go(s), and Gi{s) 
are to be expressed as 

G,(s) = ^P(k)G,(s|k) , Go(s) = ^P(k)Go(s|k), 

k k 

(7) 

and 

G{s,s')=Y,Pik)G{s,s'\k), (8) 

k 

where functions G inside the summations have the same 
meaning as the original ones under the condition of start- 
ing from a vertex of degree k. Thus, the complete solu- 
tion of the problem goes through finding the conditional 
probabilities Gi(s|k), Go(s|k) and G(s,s'|k). In the fol- 
lowing subsections we will show how to compute them in 
the general correlated case. 



A. In/Out component 

To proceed further, we first focus our attention on the 
out-component size distribution of a vertex of degree k, 
Go(s|k). Starting from a vertex of degree k = (fci, fco, kb), 
we can leave it using the ko outgoing edges and the kh 
bidirectional ones. Then, the number of reachable ver- 
tices will be the sum of the reachable vertices of each of 
the ko + kb neighbors plus 1. In mathematical terms, this 
translates into 



(9) 

I 

arrived to it from another source vertex of degree k fol- 



s s 

I 



Go(s|k) 5o(si|k) • • ■5o(sfeJk)g^(sfc„+i|k) ■ • •.g^(sfe„+feJk)(5,i+...+3,^^,^+i,^, 
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lowing one of its outgoing (bidirectional) edges. In writ- 
ing Eq.Q, we have used the fact that random networks 
are locally tree like. Equations of the type of © find in 
the discrete Laplace space their natural representation in 
terms of the generating function formalism. Using this 
formalism, these equations simplifies enormously and can 
be manipulated very easily. Within this formalism, equa- 
tion Q simplifies as 



G„(z|k) = 2[g„(z|k)]'=°[g^(z|k)] 



(10) 



where we have adopted the notation /(z) = X^s/l^)-^*- 
In what follows, we will work in the discrete Laplace 
space, using the generating function formalism. 

Functions go and satisfy the following set of coupled 
equations 



5o(z|k) = z^Po(k'lk) [.9o(^|k')]'° [g'o{^\\^')f 

k' 

gS(z|k) = z5]n(k'|k) Mz\k')f [gtiz\k')f~' 



(11) 



The term ki, — 1 in the second line of Eqs. comes 
from the fact that one of the bidirectional edges has al- 
ready been used to reach the vertex of degree k' and, 
thus, cannot be used again to leave it. Notice that this 
restriction is not needed in the first line of Eqs. Hll|) 
since, in this case, we have reached the vertex using an 
outgoing edge of the source vertex. The set of equations 
is closed for the functions g. Its solution for z = 1 
will allow us to compute the size of the GIN component, 

GIN = 1 - ^P(k) [goimf" [gtm]"' ■ (12) 

k 

The trivial solution of Eqs. (fTT|l is 5o(l|k) = 1, 
3^(1 |k) = 1, corresponding to the only case without gi- 
ant component. Therefore, the network will percolate at 
the directed level when this trivial solution becomes un- 
stable. To analyze the stability of this solution we use 
the approach adopted in and find solutions of the 
form.go(l|k) = l-e2:(k), go(l|k) = l-ey{k) in the limit 
e ^ 0. Replacing these expressions in Eqs. Hll|) and 
taking the limit e — > we obtain 



x(k) 
2/(k) 



ECk. 



x(k') 

2/(k') ; ' 



and the GIN component does not exist. In contrast, when 
A„i > 1 a non trivial solution of the set of equations 1)11(1 
exists and the GIN component emerges. 

The analysis for the in-component of individual ver- 
tices is identical to the case of the out one if we replace 
in equations (|10|) and Hll|) 



(15) 



In this case, the matrix controlling the onset of the 
GOUT component is 







9o{s\k) - 


5i(s|k) 




- 5?(s|k) 


Po(k'|k) - 


- P.(k'|k) 



*-^kk' — 



A;,'P,(k'|k) A:;P,(k'|k) 
KP,{k'\k) {ki - mik'\k) 



(16) 



As before, the condition for the appearance of the GOUT 
component is ruled by the maximum eigenvalue of the 
matrix CJ^.^., . At first glance, one could be tempted to 
conclude that, since matrices CJ^.^/ and Cj^^, are differ- 
ent, their eigenvalues are also different, leading to differ- 
ent phase transitions for the appearance of the GIN and 
GOUT components. However, it can be proved that the 
eigenvalues spectra of both matrices are identical and, 
then, both the GIN and the GOUT components appear 
simultaneously. 

It is illustrative to recover from this formalism the re- 
sults for the purely undirected and the purely directed 
cases. In indirected networks, only bidirectional edges 
are present, fc^ = and fco = 0, and the matrices turns 
into 



(17) 



recovering results in [ij]. In the case of purely directed 
networks, kb = and we obtain 



*-^kk' 

*-^kk' 



kk' 



fc;P,(k'|k) 

fc,:p,(k'|k). 



(18) 



This result generalizes the percolation theory for purely 
directed random networks developed in 0, |21| to the 



(13) case of arbitrary degree-degree correlations. 



where the matrix C^^., is defined as 



'-^kk' — 



fc;,p„(k'|k) k',Po{k'\k) 
k^Pbik'ik) {k', - i)n(k'|k) 



(14) 



The stability of the solution go(l|k) = 1, 5o(l|k) = 1 
is thus determined by the maximum eigenvalue of the 
matrix CJ^j,, , A^. When A^ < 1 this solution is stable 



B. Strongly connected component 

The analysis of the GSCC requires a more careful de- 
velopment of the ideas introduced in the previous sec- 
tion. In this case, the joint distribution of the in- and 
out-components of a vertex of degree k reads 



G{z,z'\k) = zz' [go{A'^t° [g\z,z'\k)\*'' [Uz'M 



ki 



(19) 
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where the function 17^(5, s'|k) is defined analogously to 
g'^{s\\i) and g^{s\li). It is worth to mention that, if 
the network contains bidirectional edges, G{s,s'\k) ^ 
Go{s\k)Gi{s'\li) because, in this case, such edges are com- 



mon to the in and out components of the vertex. Using 
the same reasoning as in the previous section, we can 
write down a closed equation for the joint distribution 



g\z, z'|k) = zz' n(k'|k) [go{A^)f ^'|k')] ' [ff,(z'|k')] 



-, k' -^ 



(20) 



This equation, together with Eq. H1U|I and Eqs. are 
the complete solution of the problem. The solutions for 
arbitrary values of z and z' allow to find the distribu- 
tion of the sizes of the in- and out-components of single 
vertices whereas the non-trivial solution for z = 1 and 
z' = 1, combined with Eqs. and Eqs. ifTlS)) . wiU 

provide us with a method to calculate the size of the 
different giant components of the network. 



IV. UNCORRELATED NETWORKS 

As mentioned before, real networks are usually corre- 
lated in the sense that the degrees of pairs of connected 
vertices are correlated random quantities. Nevertheless, 
uncorrelated networks are equally useful as benchmarks, 
or null models, to test topological and dynamical prop- 
erties and compare them to the results obtained in cor- 
related networks. 

When two point correlations are absent, the transition 
probabilities become independent of the degree of the 
source vertex. In this situation, after some elemental 
algebra, we obtain 



Po(k'|k) = 



, ^'^(k'lk) 



Km) 



(21) 



and 



n(k'|k) = 



ih) ■ 



(22) 



where we have made use of the fact that (kg) ~ {h). 
Using these expressions, the set of equations reduces 
to the following set of trascendent equations, 

y = {h)'^d^P{l, y,z) z = (fc,)-ia,P(l, 2/, z) 
X = {k,)-^dyP{x, 1, z') z' = {k,)-^d,,P{x, 1, z') 
z" ^{kt)-'d,„P{x,y,z"), 

(23) 

where P{x, y, z) is the generalized generating function of 
the degree distribution, that is, 

P{x,y,z)= ^ x''^y''"z'''P{h,ko,h). (24) 

The condition for the existence of a non trivial solution 
of the set of equations l|23(l is easily obtained using the 
formalism developed in the previous sections. For the 
uncorrelated case, the matrices CJ^^/ and C^^./ become 
independent of k and its maximum eigenvalue reads 



Am = 




{hko)\ , 4:{kikb)(kokb) 



{ki){h) 



(25) 



so that, whenever the condition Am > 1 is fulfilled, the 
network is in the percolated phase. As it can be seen 
from Eq. H25() , the presence of bidirectional edges alters 
the point at which the giant component arises in a non 
trivial way. The term F — (kikf,) (koki,) / (ki) (kh) mea- 
sures the strength of the coupling between directed and 
bidirectional edges. When this coupling is weak, F ~ 



and the maximum eigenvalue takes the simple form 



A„ 



max 



{kt{kb-l)) ihko) 



ih) 



(ki) 



(26) 



This result is easy to understand since, when F ^ 0, 
vertices cannot have directed and bidirectional edges si- 
multaneously and, as a consequence, the network is com- 
posed of two isolated networks, one of them purely di- 
rected and the other one containing bidirectional edges 
only. 
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In the purely undirected case, ki = and ko = and 
we recover the weh known condition for percolation in 
undirected networks with given degree distribution Eq. 

m 



{kbjkb - 1)) 



> 1. 



(27) 



In the case of purely directed networks, kt = and the 
maximum eigenvalue reads 



A. 



> 1, 



(28) 



recovering Eq. ((2Jl. 

When r 3> 1, Am must be computed using Ea. l(^ 
and, in general, will depend on the density of edges, as 
well as on the type of correlations between directed and 
bidirectional edges. In particular, a positive correlation 
between kb and ki or kg can strongly favor the emergence 
of the giant component even if the density of bidirectional 
edges is very small. We will illustrate this point in the 
example of the next section. 

To finish the uncorrelated analysis, let us com- 
pute the relative sizes of the giant components. Let 
{xcVci Zc, z'^, z'^) be the non trivial solution of the set 
ofEqs. (123) , then using Eq. ((SJ and Eq. ©, the relative 
sizes of the different giant components of the network 
read 

GIN - 1 - P(l, y„ zc) , GOUT = 1 - F(xe, 1, 4) (29) 



GSGC = 1 - F(l, yc, z,) - P{x„ 1, 4) + P{x„ y„ z'^) 

(30) 

In the next section, we present a practical application 
of this formalism. 



V. BIDIRECTIONAL EDGES AS 
PERCOLATION CATALYSTS 

Suppose we have a purely directed network with de- 
gree distribution P{ki, ko) and no two point degree cor- 
relations. Suppose also that the network is in a regime 
in which the GWCC exists but not the GSCC, that is. 



{ki ko) 



< 1 and 



{{h + ko){k^ + ko - 1)) 

{{h + ko)) 



> 1. 



(31) 



Now we transform the original network by converting 
each directed edge into a bidirectional one with probabil- 
ity p. After this transformation, we end up with a net- 
work with pE bidirectional edges and (1 — p)E directed 
ones, where E is the original number of directed edges on 
the network. The degree distribution of the transformed 
network can be written, in the discrete Laplace space, as 

Pp(zj,Zo, Zb) = P{pzb + (1 -p)z^,pz + (1 -p)zo). (32) 




FIG. 2: Relative size of the giant strongly connected compo- 
nent as a function of the conversion probability p. Simulation 
results are for a single network with A'^ — 10^ vertices. Solid 
lines are the numerical solution of the set of equations 12311 . 




FIG. 3: Relative size of the giant in component as a function 
of the conversion probability p. Simulation results are for a 
single network with TV = lO" vertices. Solid lines are the 
numerical solution of the set of equations I23II . 



This transformation undergoes a phase transition as we 
increase the value of p. When p = 0, the network is 
purely directed and, by construction, it has no GSCC. 
Whenp = 1, all edges become bidirectional and, thus, the 
GSCC is identical to the GWCC of the original network. 
Therefore, at some intermediate value Pc, the network 
percolates and a GSCC emerges. The value of Pc can 
be easily obtained using the expression for the maximum 
eigenvalue, Eq. (|25|) . and the final degree distribution, 
Eq.dSa). 

The most interesting case corresponds to networks 
with marginal degree distributions following power laws 
of the form Pi{ki) ~ k~'^^ and Po{ko) ~ ko^" with 
li,lo < 3. When the transformation described above 
is performed in this type of networks, some of the terms 
in Ea. H25|) are proportional to {kf) and (fc^) and, con- 
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sequently, Km — > oo in the thermodynamic hmit. This, 
in turn, imphes that Pc — 0, that is, even an infinitely 
small fraction of bidirectional edges suffices to percolate 
the network. 



A. Numerical simulation 

The check the accuracy of our theory in the case of 
power law marginal degree distributions, we have per- 
formed extensive numerical simulations. We first gener- 
ate purely directed random networks with degree distri- 
bution of the form P{ki,ko) — Pi{ki)Po{ko) and no two 
point correlations. The in and out degree distributions 
are taken to be identical and to follow a scale-free form 
of the type 

P,{k) = Po{k)= \ (33) 



where C(7) is the Zeta Riemann function. To generate 
a purely directed random network, we use a natural ex- 
tension of the configuration model 0, 25, 26] -an al- 
gorithm intended to generate uncorrelated random net- 
works with a given degree distribution. The algorithm 
starts by first assigning to a set of N vertices a pair of 
"stubs" , one of them incoming and the other one out- 
going, ki and /cq, randomly drawn from the distribution 
P{ki, ko). The only requirement is that ki — fco, 
whenever one wishes to close the network. The network is 
constructed by selecting pairs of in and out stubs chosen 
uniformly at random to create directed edges, avoiding 
multiple, bidirected, and self-connections among vertices. 
Once the network has been assembled, each directed edge 
is transformed into a bidirectional one with probability 
p, and the relative sizes of the giant components are mea- 
sured. 

Figures |21 and O show simulation results for a scale- 
free network following Ea. H33l) with exponent 7 = 3 and 
size N = 10® as compared to the numerical solution of 
Eas. H23|l . As it can be seen, the agreement between sim- 
ulation results and the theoretical prediction is excellent. 
The relative sizes of the GSCC and the GIN are shown 
as a function of the conversion probability p for different 
values of the distribution parameter Pqj the probability 
of nodes having null in or out degree. Even for very small 



values of p, the GSCC and the GIN are evident. As ex- 
pected, small values of Pq favor the growth of bigger giant 
components. Then, bidirectional edges act as a percola- 
tion catalysts, favoring the appearance of a fine structure 
in the giant connected component. The scale- free prop- 
erty of many real networks is, once more, indicative of 
interesting features, since, in this case, the presence of 
an infinitesimal fraction of bidirectional edges is enough 
to ensure percolation at the level of the directed compo- 
nents . 

VI. CONCLUSIONS 

We have derived a very general formulation of the 
theory of percolation in directed random networks with 
bidirectional edges and arbitrary two point degree cor- 
relations. Our formalism accounts for all the previ- 
ously known results for percolation in purely directed 
and purely undirected random networks, which stand as 
limiting cases of our theory. The percolation threshold 
for the most general situation is derived as a function of 
the maximum eigenvalue of the connectivity matrices. In 
particular, for networks with no two point correlations, 
explicit expressions are provided depending on the first 
and second moments of the degree distribution P(k). In 
this case, we have also shown that bidirectional edges act 
as a catalyst for percolation, favoring the emergence of 
the GSCC, and for scale- free networks, only an infinites- 
imal fraction of bidirectional edges is needed. 

After the completion of this work, we have become 
aware of a recent preprint '27| where the classical 
Susceptible-Infected-Recovered (SIR) model of epidemi- 
ology is analyzed in uncorrelated directed random net- 
works with bidirectional edges. Since there exist a map- 
ping between the SIR model and percolation theory, some 
of the results derived in that reference overlap our results 
for the uncorrelated case. 
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